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Abstract 

A relativistic framework for computing the nuclear transparency extracted from 
A{p, 2p) scattering processes is presented. The model accounts for the initial- and 
final-state interactions (IFSI) within the relativistic multiple-scattering Glauber ap- 
proximation (RMSGA). For the description of color transparency, two existing mod- 
els are used. The nuclear filtering mechanism is implemented as a possible explana- 
tion for the oscillatory energy dependence of the transparency. Results are presented 
for the target nuclei 7 Li, 12 C, 27 A1, and 63 Cu. An approximated, computationally 
less intensive version of the RMSGA framework is found to be sufficiently accurate 
for the calculation of the nuclear transparency. After including the nuclear filtering 
and color transparency mechanisms, our calculations are in acceptable agreement 
with the data. 
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1 Introduction 



The transition region between nucleon-meson (hadronic) and quark-gluon 
(partonic) degrees of freedom is a topic of longstanding interest in nuclear 
physics. A promising observable to map this transition is the transparency of 
the nuclear medium to the propagation of hadrons. In A(p, 2p) experiments, 
the nuclear transparency is defined as the ratio of the cross section per nu- 
cleon to the hydrogen one. Accordingly, the nuclear transparency is a measure 
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for the attenuation effects of the spectator nucleons on the impinging and 
outgoing protons. 

In the conventional Glauber picture [Tj, the nuclear transparency extracted 
from A(p, 2p) reactions is predicted to be rather constant for incoming mo- 
mentum larger than a few GeV/c. The color transparency (CT) phenomenon 
suggests an anomalously large transmission probability of protons through 
nuclei [2113] and leads to a nuclear transparency that increases with incoming 
momentum. The experimental A(p, 2p) data [4.5,6] suggest a CT-like increase 
in the transparency for impinging proton momenta between 5 and 10 GeV/c. 
For higher momenta the nuclear transparency falls back to the Glauber level. 
This oscillatory energy dependence is not unique to the A(p, 2p) nuclear trans- 
parency: it has been observed or hinted at in pp elastic scattering [7], elastic 
Tip fixed-angle scattering [8f5] . pion photoproduction PfTU] . and deuteron pho- 
todisintegration [TTlfTS] . 

One possible interpretation of this energy dependence of the A(p, 2p) trans- 
parency is provided by the presence of two terms in the free pp scattering 
amplitude [T5ffT4] . Ralston and Pire [T3] suggested a combination of the fol- 
lowing two components. First, the quark-counting component which follows 
the dimensional scaling law [15] and represents a small object (a point-like 
configuration, PLC). Second, the Landshoff component [16] which is asso- 
ciated with normal-sized configurations. The interference between these two 
amplitudes induces the oscillation of the pp cross section about the scaling 
behavior. Inside the nuclear medium, the Landshoff component will be sup- 
pressed because of the strong interactions with the spectator nucleons, while 
the quark-counting component will escape the nucleus with relatively small 
attenuation due to CT. This phenomenon is called "nuclear filtering" (NF) 
[T7] : the nucleus filters away the normal-sized components in the hadron wave 
functions. Accordingly, the nucleus plays an active role in selecting small-sized 
components. The NF model is able to reproduce qualitatively the observed 
bump in the A(p, 2p) transparency. 

Recently, a relativistic and cross-section factorized framework for A(p, 2p) re- 
actions has been proposed [18] . In this Letter, this formalism is extended to 
incorporate the Ralston-Pire model for the pp scattering amplitude and the 
concepts of CT and NF. For the description of the initial- and final-state in- 
teractions (IFSI) of the impinging and two outgoing protons, we adopt the 
relativistic multiple-scattering Glauber approximation (RMSGA) [18] . This 
relativistic extension of the Glauber model, which was originally developed 
to describe A(e,e'p) observables [TSPU] . models the IFSI as consecutive cu- 
mulative scatterings with the individual spectator nucleons in the nucleus. 
To compute the effects of CT, we consider the quantum diffusion model of 
Ref. [21] and an alternative treatment presented in Ref. [22]. The comparison 
between the two CT models is made in a consistent way. 
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This Letter is organized as follows. In Section [2j the RMSGA A(p, 2p) formal- 
ism of Ref. [18] is extended to take into account both the quark-counting and 
the Landshoff contribution to the hard pp scattering amplitude. A numerically 
convenient approximate form of the RMSGA framework is introduced. Fur- 
ther, we describe two different methods to deal with CT, thereby indicating 
the similarities and differences. Section [3] presents the nuclear transparency 
results for the target nuclei 7 Li, 12 C, 27 Al, and 63 Cu. The effect of IFSI, CT 
and NF is discussed. Furthermore, the accuracy of the approximated RMSGA 
approach is investigated and the two CT treatments are compared. Finally, 
Section @] summarizes our findings and states our conclusions. 



2 Formalism 



In the Ralston-Pire approach [23] , the spin-averaged pp scattering amplitude 
consists of the quark- counting (QC) and the Landshoff (L) contribution: 

M» = M% + M?. (1) 



The Ai q C term results in a differential cross section that scales like da pp / dt oc 
s~ 10 and the Landshoff term can be related to it through 

M " = f VlG^ er*(*M+«0 M% c . (2) 



Here, s and t are the Mandelstam variables. Further, p\ = 0.08, <ft i s ) = 
In | In 01 Q eV 2 }, and 5\ = —2.0. These values were determined from a 
fit to the pp data at 90° [23] . The sole parameter which remains undetermined 
is the sign of the phase difference (s) + 5i between the quark- counting and 
the Landshoff term. Therefore, both signs will be used in the calculations. 

Incorporating the Ralston-Pire approach and the NF mechanism into the 
A(p, 2p) formalism of Ref. [15] . the amplitude for the p (E p i,pi, m s u)+A (e a , k^, 0" 

p [E kl , ki, m s if) + p [E k2 , k 2 , m s2f ) + A-l (E a _ x , k A _ u J R M R j reaction be- 
comes 



Mff p) = £ (M% c ) u(p m , m s )0j MSG ^ CT {p m ) 

m 3 m s ii,m s ,m sl f,m s2 f 1 

+ E (Ml P ) mBli , ms ,m 8lf ,m s2f "(Pm, m s )^ SGA (p m ) , (3) 

where p m = k\ + k 2 — p± is the missing momentum and a,\ refers to the state 
wherein the struck proton resided. In this expression, the effect of IFSI is ac- 
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counted for through the distorted momentum-space wave functions 0R^sga+ct 
and 0^ MSGA (p m ). Since the quark-counting term is associated with PLCs, the 
corresponding momentum-space wave function 0Rmsga+ct includes the 
effect of CT. The Landshoff term, on the other hand, corresponds with a 
hadron of normal size. Consequently, the IFSI can be computed in standard 
Glauber theory. Using the spin-averaged pp matrix element of Eq. (JT]), the 
squared A(p, 2p) matrix element for knockout from the ot\ shell can be cast in 
the form 



if 



(P,2p) 



RMSGA+CT 



(Pn 



+2Re 



M^ c (Miru(p m ,m s ) c 



RMSGA+CT 



(Pr, 



(u(p ri 



+ \M P L P \ 2 u{p m ,m s )<P^{p m ) 



iRMSGA 



RMSGA 



iPr, 



(4) 



with m the struck nucleon's generalized angular momentum quantum number. 
The differential cross section is obtained as an incoherent sum of the squared 
matrix elements over all proton levels «i, thereby factoring in the occupation 
number of every level. The momentum-space wave function is defined as [18] 



RMSGA( + CT) ^ 



(5) 



where the relativistic bound-state wave function <p ai (r) is computed in the 
Hartree approximation to the a — uo model |24j, using the Wl parametrization 
for the different field strengths [25]. The t 5 I ^ [ I SGA( ' +CT ' ) (r) operator accounts 
for the IFSI effects and is the subject of the forthcoming discussion. Hereafter, 
results obtained on the basis of Eq. (0} are dubbed RMSGA+CT+NF. 

In the RMSGA framework, the IFSI operator takes on the form [18] 




where r denotes the point of collision between the struck and incoming pro- 
ton, and fj (j = 2, . . . , A) are the positions of the frozen spectator nucleons. 
Further, the z axis lies along pi, z' along and z" along k 2 - The b, b', and b" 
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planes are perpendicular to these proton momenta. Reflecting the diffractive 
nature of pN collisions at GeV energies, the profile function for pN scattering 
is parametrized as 



, /, = qffi 0) (1 - ie pN (AQ) | u- | 



The total pN cross section cr*^ (k), the slope parameter /3 p at (k), and the ratio 
of the real to the imaginary part of the scattering amplitude e P N (k) depend 
on the proton momentum k. In the numerical calculations, their values are ob- 
tained through interpolation of the iViV scattering database from the Particle 
Data Group [26J. 

Eq. ([6]) is a genuine A-body operator and the integration over the coordinates 
of the spectator nucleons makes its numerical evaluation very challenging. 
The IFSI operator can be rewritten in a numerically more convenient form by 
adopting the following approximations. First, the squared wave functions of 
the spectator protons (neutrons) are replaced by 1/ [Z — 1) [1/N) times the 
proton (neutron) density of the residual nucleus. Second, one assumes that 
these densities are slowly varying functions of b, while T p n (k, bj is sharply 

peaked at b = 0. This allows one to approximate the IFSI operator of Eq. <^ 
by the one-body operator [27J: 



£RMSGA' ^ _ -Q e -50-*$Oi)(l-iepJ\r(Pi)) Jl^^j PN&Zj) 



N=p,n 

x e 



X e -5<w(^)(l-V(fe)} XT<i^"w(6",%") _ 

Here, p p and p n are the proton and neutron density of the residual nucleus. 
They reflect the spatial distribution of the scattering centers inside this nu- 
cleus. Henceforth, calculations based on Eq. (jSJ) are labeled as RMSGA'. 

The essential assumption of CT is that the impinging proton compresses to a 
PLC as it hits a target nucleon, after which the outgoing protons expand from 
PLCs to normal-sized objects as they move through the nucleus. To account 
for the reduced interaction of a PLC with the nuclear medium, the total cross 
sections in Eqs. ((71) and (jHD are replaced by effective ones. In the partonic 
model of Farrar et al. [21] (denoted by FLFS), the interaction cross section is 
argued to be 



,21,2 

FLFS 

pN 



(P, Z) = a^{[f- + ^(l - ^ )}9(l h - Z) + 6(Z -l h )). (9) 
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Here, Z is the distance from the hard interaction point along the trajectory of 
the particle, n = 3 is the number of constituents in the proton, and (k 2 ) 1 ^ 2 = 
0.35 GeV/c is the average transverse momentum of a parton in a hadron. The 
quantity 1% ~ 2p/AM 2 is the hadronic expansion length, i.e., the propagation 
distance at which an expanding hadron reaches its normal hadronic size, and 
depends on the hadron momentum p and the squared mass difference AM 2 
between the intermediate PLC and the normal-sized hadron. It is commonly 
assumed that 0.7 < AM 2 < 1.1 (GeV/c 2 ) 2 are reasonable values. 

Starting from a hadronic picture, Jennings and Miller [22] suggested the fol- 
lowing alternative expression for the effective cross section 



^(p,Z) = alt(l-^e^ z ), (10) 



with p the proton momentum and p* the momentum of a baryon resonance 
with a complex mass M* and the same energy as the nucleon, i.e., (p*) 2 = 
p 2 + M 2 — (M*) 2 . The expression for the effective cross section emanates 
from the intermediate PLC being a superposition of the nucleon ground state 
and a nucleon resonance. The imaginary part of M* ensures the decay of the 
intermediate state to an asymptotically free, normal-sized proton. Like the 
FLFS approach of Eq. (jSJ), Eq. ( fTOl) considers one excited state in the PLC. 
It is worth noting that both the FLFS and JM model take into account the 
suppression of interaction in the collision point and the time evolution of the 
PLC to a normal-sized proton during its propagation through the nucleus. 

In our numerical calculations, we will also consider the standard RMSGA+CT 
picture. In this scenario, the entire wave packet of the incoming and outgo- 
ing protons is assumed to be in a PLC, which propagates through a passive 
nuclear medium. This amounts to neglecting the Landshoff term in the ampli- 
tudes M pp and Mff p) of Eqs. © and ©. Finally, in the standard RMSGA 
calculations, both the Landshoff term and CT effects are neglected. 

For the free pp scattering cross section da pp /dt, the parametrization as pre- 
sented in Ref. [28] is used. This parametrization combines the 6 c .m. dependence 
suggested by [29] and the Ralston-Pire separation of Eq. (T5]). 
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3 Nuclear transparency results 



The nuclear transparency is computed as the ratio of the cross sections in- 
cluding and excluding IFSI effects: 

°RPWIA 



The relativistic plane wave approximation (RPWIA) limit is reached by setting 
the IFSI operator 5ifsi if) equal to one in Eq. and corresponds with a 
calculation which ignores IFSI. The numerator and denominator of Eq. (TTTT) 
are obtained by integrating the corresponding differential cross sections over 
the phase space defined by the kinematic cuts. In our calculations, we adopted 
identical cuts as in the experiments jlfofo] and assumed a flat experimental 
acceptance within the kinematical ranges for each data point. 

We consider the experimental nuclear transparency values as presented in |30j . 
The incident lab momentum varies from 5.9 to 14.4 GeV/c and the scattering 
angle is near 90° in the pp center of mass. The Mandelstam variable \t\ ~ 
(s - 4M 2 ) /2 extends from 4.7 to 12.7 (GeV/c) 2 . In Ref. [30], the originally 
published values of Carroll et al. [I] were rescaled using the improved nuclear 
momentum distributions of Ref. [21], thereby making them consistent with 
the data of Refs. [5llo] . 

First, we address the energy dependence of the l2 C(p, 2p) transparency and 
study the role of IFSI, CT and NF. Further, we use the 12 C(p, 2p) calculations 
as a test case to determine the accuracy of the IFSI operator of Eq. (jSJ) rel- 
ative to the expression of Eq. (jBj). Fig. [TJ displays the 12 C transparency as a 
function of the incoming proton momentum p\. The solid curves represent the 
full RMSGA calculations, whereas the RMSGA' results are shown as dashed 
curves. Three different scenarios were considered. As expected, the standard 
RMSGA calculations lead to a nuclear transparency that is almost indepen- 
dent of the beam momentum. The main effect of the IFSI is to reduce the 
nuclear transparency from the asymptotic value of 1 to ~ 0.15. The inclu- 
sion of CT effects produces a transparency linearly rising with energy. The 
increase relative to the RMSGA result is highly dependent on the adopted 
model and corresponding parameters for CT. The curves including CT shown 
in Fig. □ adopt the FLFS model with AM 2 = 0.7 (GeV/c 2 ) 2 . The increase of 
the transparency is consistent with the data in the range 5-10 GeV/c, but the 
RMSGA+CT picture fails to explain the drop in the transparency at higher 
momenta. Our RMSGA and RMSGA+CT predictions confirm the results of 
[32] . A better agreement with the data is obtained when adding the mechanism 
of NF. Compared to the RMSGA+CT results, the transparency is increased 
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at intermediate momenta (5-10 GeV/c) and decreased at higher momenta, 
two effects which improve the description of the data. A similar result was 
obtained in Ref. [33] where the JM model of CT was used. 

Concerning the comparison of the RMSGA and RMSGA' results, it can be 
inferred from Fig. [I] that both approaches yield nearly identical results which 
differ at the 2-3% level. Consequently, the operator of Eq. (jSJ) is considered 
sufficiently accurate for the calculation of the nuclear transparency and will be 
used in the remainder of this work. We wish to stress that the computational 
cost of Eq. (jSJ) is about a factor of 10 3 lower than the full-blown RMSGA 
operator of Eq. (J6j). 

Figs. [2] and [3] are devoted to a comparison of the different CT models. Re- 
sults of the FLFS quantum diffusion model are plotted for AM 2 = 0.7 and 
1.1 (GeV/c 2 ) 2 . For the M* parameter of the JM model we consider three dif- 
ferent values, representing the A, the Roper resonance, and the average of 
the lowest P-wave N* resonances. For the imaginary part of M* a value of 
150 MeV was taken. The lowest values of the parameters AM 2 and M* in- 
duce the strongest increase of the RMSGA+CT transparency with the beam 
momentum pi and also lead to the largest deviations between the predictions 
including the NF mechanism and the corresponding RMSGA+CT results. 

Fig. [2] shows that after including NF and CT, the calculations correctly re- 
produce the maximum in the 12 C(p,2p) transparency at about 9.5 GeV/c, 
but badly fail to fall back low enough to account for the 14.4 GeV/c data 
point. Ralston and Pire [13J, on the other hand, do succeed in reproducing the 
maxima and minima in the nuclear transparency data. However, they assume 
that the transparency of the quark-counting term is beam-energy indepen- 
dent, a rather speculative assumption. The CT-induced increase of the quark- 
counting transparency with energy causes our RMSGA+CT+NF predictions 
to rise again at a momentum pi ~ 12 GeV/c. 

In the FLFS as well as the JM approach, the RMSGA+CT+NF predictions 
reproduce the general trend of the data, but no variant achieves very good 
agreement. Furthermore, it is not possible to unambiguously determine the 
value of the parameters AM 2 and M*, as the best choice for these parameters 
depends on the target nucleus under consideration. Fig. [2] suggests that for 
the 12 C(p, 2p) reaction M* = (1440 + 150z) MeV leads to the best agreement, 
while for the AM 2 parameter no "best" choice can be put forward. Fig. [31 
on the other hand, shows that for the 27 A1 target nucleus the FLFS results 
are systematically below the data in the region below 10 GeV/c. Using M* = 
(1440 + 150i) MeV in the JM model also does not increase the transparency 
high enough so as to match the 6 and 10 GeV/c data points, only with M* = 
(1232 + 150i) MeV is the CT-induced increase of the transparency strong 
enough. 
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None of the results shown in this Letter include the color screening effect (CSE) 
[34] . This QCD effect suggests the suppression of the small-size configurations 
in bound nucleons. We deem that the CT parameters are so badly constrained 
that controlling additional mechanisms is out of reach for the moment. In both 
the RMSGA+CT and the RMSGA+CT+NF calculations, the inclusion of the 
CSE decreases the transparency by 6-12%, with the largest effect occuring at 
higher p\. 

Another effect that can be studied in Figs. [2] and [3] is the influence of the 
sign of <p (s) + 5i on the RMSGA+CT+NF results. For the FLFS model of 
CT, the differences between calculations using both signs of (s) + 8\ are 
minor. As already observed by Jennings and Miller [33], the results using the 
JM model of CT are rather sensitive to this sign. The discrepancy between 
the FLFS- and JM-based calculations arises from the different structure of 
the effective cross sections ([HD and (fTDl . Indeed, these effective cross sections 
not only determine the attenuation of the quark-counting term in the nuclear 
medium, but also the phase difference between the quark-counting and the 
Landshoff term. Whereas the real parts of both effective cross sections are 
quite similar, the FLFS effective cross section (jH]) is purely real, while its 
JM counterpart also has an imaginary part. This imaginary part causes the 
enhanced sensitivity of the JM results to the sign of <fi (s)+5\. As for which sign 
causes the best agreement with the data, no firm conclusions can be drawn. 
Indeed, the optimum choice for the sign of (s) + 5\ depends on the used 
CT model (FLFS or JM), the value of the parameter AM 2 or M*, and the 
target nucleus. For the 7 Li(p, 2p) reaction, the positive sign provides the best 
agreement, while the 63 Cu transparency data rather require a negative sign. 

The A dependence of the nuclear transparency at two values of the incoming 
momentum pi is studied in Fig. HI The standard RMSGA calculations fall 
considerably below the data. Further, none of the RMSGA+CT+NF calcula- 
tions succeed in simultaneously describing the data for all target nuclei. While 
the FLFS approach agrees with the 7 Li and 12 C data points rather well using 
AM 2 = 0.7 (GeV/c 2 ) 2 , the FLFS results tend to underestimate the 27 Al and 
63 Cu data. With regard to the M* parameter of the JM model, for the 7 Li 
and 12 C nuclei a value of M* = (1440 + 150*) MeV seems acceptable, whereas 
the heavier 27 Al and 63 Cu nuclei need a smaller M* value. A general feature 
of the RMSGA+CT+NF predictions is that their A dependence is steeper 
than the data. Finally, the dot-dashed and dotted curves indicate that the 
6 GeV/c data are proportional to A~ 2 / 3 , while at 10 GeV/c the A dependence 
of the data is more gradual (T oc A~ l l 3 ). This trend is not reproduced by 
the standard RMSGA predictions, which are almost independent of the in- 
coming momentum. When CT and NF effects are included, the softening of 
the A dependence with increasing incoming momentum is also present in the 
calculations. For example, the FLFS calculations with AM 2 = 0.7 (GeV/c 2 ) 2 
correspond to T oc A~ om at 6 GeV/c incoming momentum and T oc A~ ' 70 
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at 10 GeV/c. 



4 Conclusions 

In conclusion, we have developed a relativistic framework to calculate the 
nuclear transparency for A(p, 2p) processes. A relativistic multiple-scattering 
Glauber model was used to account for the IFSI. To reduce the computational 
cost of the RMSGA calculations, some additional approximations were made. 
The predictions with the full and the approximated RMSGA approach agree 
at the few percent level. Thus, to determine integrated quantities such as 
the nuclear transparency, a valid alternative for the computationally intensive 
RMSGA framework is available. 

Using the concept of nuclear filtering, our calculations are in qualitative agree- 
ment with the data, thereby confirming earlier calculations [33 35 36J. Further- 
more, our calculations seem to indicate that CT is imperative to increase the 
calculations to the level of the data. The same conclusion was reported in 
Refs. [3"oTI3"T] . The quantitative description of the data, however, is far from 
perfect and it is not possible to constrain the magnitude of the parameters in 
the CT models. 

The indications for CT in A(p, 2p) reactions are not necessarily in contra- 
diction with the results from A(e, e'p) experiments. Although the A(e, e'p) 
transparencies show no significant increase with the four-momentum trans- 
fer Q 2 [38] and can be reasonably reproduced in the RMSGA framework, 
the existence of CT can not be excluded since the predicted effect is small 
|39j . This is caused by the small expansion times of the PLC to normal size 
at the present A(e, e'p) kinematics. The effect of CT is more pronounced in 
the A(p, 2p) transparency for different reasons. First, in the A(p, 2p) reac- 
tion there are three particles that can experience CT instead of only one in 
A(e, e'p) reactions. Second, A(p, 2p) data are available up to Q 2 = \t\ values 
of 12.7 (GeV/c) 2 , while A(e,e'p) transparency experiments are restricted to 
Q 2 < 8 (GeV/c) 2 . 

A number of uncertainties involving the A(p, 2p) transparency remain. One 
subject of discussion is the size of the Landshoff term. According to Botts 
et al. jlD], this term might also be small-sized. This would make the sur- 
vival probability of the Landshoff and the quark-counting term rather similar 
and would weaken the oscillations in the energy dependence of the computed 
transparencies. Apart from the Ralston-Pire picture discussed above, other ex- 
planations of the energy dependence of the transparency have been suggested. 
For example, Brodsky and de Teramond [H] interpreted the oscillatory be- 
havior in terms of two broad baryon resonances associated with strange and 
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charmed particle production thresholds, interfering with a perturbative QCD 
background. An improved set of data, particularly at higher energies, is es- 
sential to clarify these issues. The 50 GeV proton synchrotron that is under 
construction at J-PARC [JT] opens great opportunities for this research. 
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Fig. 1. The nuclear transparency for the 12 C(p,2p) reaction as a function of 
the incoming lab momentum p\. The full RMSGA (solid lines) are compared to 
the RMSGA' (dashed lines) results. The different curves represent the RMSGA, 
RMSGA+CT and RMSGA+CT+NF calculations. The CT effects are calculated 
in the FLFS model [21] with AM 2 = 0.7 (GeV/c 2 ) 2 and the results including the 
mechanism of NF are obtained using the positive sign of 4> (s) + Si. Data are from 
Refs. |4T5T6] . 
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Fig. 2. The 12 C(p, 2p) transparency versus the incoming lab momentum p\. The 
upper (lower) panel depicts results using the FLFS (JM) model for CT. Calculations 
including the effects of CT and NF with the positive (solid lines) and negative 
(dashed lines) sign for <fi (s) + 5i are shown, along with the RMSGA+CT predictions 
(dot-dashed lines). Data are from Refs. [I] (solid triangles), [5] (solid squares), and 
[6] (solid circles). 
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Fig. 3. As in Fig. [U but for 27 A1. Data are from Ref. @]. 
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Fig. 4. A dependence of the nuclear transparency at two values of the incoming 
lab momentum p\. The standard RMSGA calculations are represented by dashed 
curves, while the solid curves are RMSGA+CT+NF calculations with the positive 
sign for 4>(s) + 6\. The solid curves correspond with different descriptions of the 
CT effects, as indicated by the legend. The dot-dashed (dotted) curves display the 
A^ 1 / 3 (A~ 2 / 3 ) parametrization, normalized to the 7 Li data points. Data are from 
Ref. @]. 
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